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Abstract 

We provide a necessary and sufficient condition for tlie ratio of two 
jointly Q-Frecliet random variables to be regularly varying. This condition 
is based on the spectral representation of the joint distribution and is 
easy to check in practice. Our result motivates the notion of the ratio tail 
index, which quantifies dependence features that are not characterized by 
the tail dependence index. As an application, we derive t he asymptotic 
behavior of the quotient correlation coefficient proposed in lZhan^ (|2008l ) 
in the dependent case. Our result also serves as an example of a new 
ty pe of regular variatio n of products, different from the ones investigated 
bv lMaulik et"ai] l|2002l ). 

Keywords: Regular variation, multivariate Frechet distribution, spectral 
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1 Introduction 

Regular variation is often used to describe the tail behavior of random variables. 
A measurable function U : IR+ — > M-|_ is regularly varying at infinity with index 
p, written U G RV-p, if 

U{tx) 

lim — -— — X , for all a; > . 

t^oo U{t) 

When p = we say U is slowly varying. A random variable Z is regularly varying 
with index p > 0, written Z e RV-p, if th e tail dis t ributi on = P(Z > t) 
is regularly varying: F S RV-p (see e.g. iResnidd (119871 )). For the sake of 



simplicity, we only consider non-negative random variables. 

In extreme value theory, the notion of regular variation plays an important 
role in characterizing the domain of attraction of random variables. Namely, a 
random variable Z is regularly varying with index a > 0, if and only if 
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as n — >■ cxD, where {KnjnGN is some normalizing sequence, Zi, ^2, . . . are indepen- 
dent and identically distributed (i.i.d.) copies of stands for convergence 
in distribution and Ca is a standard a-Frechet random variable with distribution 

= P(Ca < = exp(-i-") , i > . 

See e.g. iResnick ()l987l) . P roposition 1.11. In this case, the random variable Z 
is said to be in the domain of attraction of $q,. The notion of regular variation 
for multivariate random ve ctors and stochastic p rocesses have also been i nvest i- 
gated extensively. See e. g. Resnick ( 1987 . 2007), jie Haa ii and Ferreira ( 200d ). 
Balkema and Embrechts (2007), and Hult and Lindsko3 (|2005l ). among others. 
Other applications of regular variation include, just to menti on a few , the d o- 
main of attraction of partial sums of i.i.d. random v ectors (Rvaceva (|l962f) ). 
large deviations of regularly varying random walks (|Hult et al.l (|2005[ )) and 
finite-dimension al distribution s of t he stat i onary solution of stochastic recur- 
rence equations (jKestenl (Il973l) and lGoldid (|l99l[ )l. 

In this paper, we consider the regular variation of ratios of two random 
variables X and Y . From a practical point of view, the ratio XjY can be seen 
as a random normalization of X by Y . In extreme value theory, when modeling 
extremal behaviors, certain normalizations (thresholding) of values are often 
required. Random normalization sometimes has appealing theoretical propertie s 
and simplifies the statistical applications (see e.g. lHeffernan and Resnickl (2007), 
Section 4). 

By viewing XjY as the product of X and 1/F, the problem is closely related 
to the regular variation of products of random variables. Whe n the two random 
variables are indepe ndent, this problem has been addressed bv iBreimanI ( 19651 ). 
On the other hand, iMaulik et al.l (|2002i ) investigated certain dependence cases, 
which were then applied to the modeling of network traffic. 

We address the regular variation of ratios in a specific case. Namely, we sup- 
pose that (X^Y) has a bivariate a-Frechet distribution (or, {X,Y) are jointly 
a-Frechet), i.e., for all a,b > 0, max{aX,bY) has a-Frechet distribution. To 
study the bivariate g-Frec het distribution s, an efficient tool i s the spectral rep- 
resent ation introduced bv Ide Haan ( 19841) and developed bv IStoev and Taggu 
( 2005h (a brief review will be given in Section^. Based on the spectral represen- 
tation, we provide a necessary and sufficient condition for X/Y to be regularly 
varying (Theorem [1} . If this is the case, then the regular variation index is 
referred to as the ratio tail index of X/Y. We demonstrate that our condition 
is easy to check through a few popular models. 

Our specific se ttin g provides examples that are not covered by the results 



BreimanI ( 19651 ) and iMaulik et al 



(|2002l ). Furthermore, we show that the ra- 
tio tail index does not characterize the dependence between X and Y in the tra- 
ditional sen se. We will compare the ratio tail index and the tail dependence in- 
dex (s ee e.g. Sibuva ( 1960l ). de Haan and Resnick (1977) and Ledford and Tawn 
( 1999) ). which has been widely used to quantify asymptotic (in) dependence of 
random variables. 

As the main application of our result, we derive the asymptotic behavior 
of the quotient correlation coefficient (Theorem [2|) for jointly Frechet distribu- 
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tions. This coefficient was proposed by IZhand (l2008h and applied to test the 
independence of two random variables X and Y . In this so-called gamma test, 
the quotient correlation coefficient is based on the independent samples of ra- 
tios X/Y and Y/X. The asymptotic behavior of the coefficient has so far been 
studied only in the case when X and Y are independent Frechet. 

Our result provides new theoretical support for the gamma test. We show 
that, when {X, Y) is bivariate Frechet, the power of the gamma test is high most 
of the time. Indeed, the asymptotic behavior of the quotient correlation coeffi- 
cient is essentially determined by the ratio tail indices of {X,Y) and (Y,X), if 
they exist. Furthermore, if the ratios have lighter tails than the single variables, 
then the gamma test rejects the null hypothesis with probability going to one as 
the sample size increases to infinity (Corollary [3]). We also show that, when the 
ratios have tails equivalent to the ones of marginals, then in a 'worst scenario', 
the gamma test performs also reasonably well (Example [5]) . 

The paper is organized as follows. Section [5] provides preliminaries on the 
regular variation and the spectral representations of bivariate a- Frechet random 
vectors. Section[3]is devoted to the characterization of regular variation of X/Y, 
based on the joint distribution of {X/Y,Y). Section |4] calculates the ratio tail 
indices for several well-known examples. In Section[Sl we review Zhang's gamma 
test and apply the result in Section [3] to derive the asymptotic behavior of the 
quotient correlation coefficient in the dependent case. Sectio n. [5] provides a proo f 
of the joint distribution of {X/Y,Y), based on a result by IWeintraubl (Il99ll) . 
Finally, in Section [7] we briefly discuss the connection between our results and 
some related works. 



2 Preliminaries 

In this section, we review the spectral representation of bivariate a-Frechet 
distributions. We also introduce some notations that will be used in the rest of 
the paper. 

We will focus on bivariate 1-Frechet random vector {X, Y). Every bivariate 
1-Frechet random vector {X, Y) has the following spectral representation: 

{X,Y) ^ (7 f{s)M{dslf gis)Mids)) . (1) 

Here, ^ stands for the extremal integral, {S,Bs,lJ^) is a standard Lebesgue 
space, for example, a Polish space (S, p) with a cr-finite measure /i on its Borel 
sets Bs, f,g are measurable non- negative and integrable functions on (S, Bs, n), 
and M is a 1-Freche t ran dom sup-measure on (S. B s) with control measure fj, (see 
e.g. de HaanI (|l984l ) and lStoev and TaaaiJ ^QM))- The functions / and g are 



called the spectral functions of X and Y, respectively, and the joint distribution 
of {X, Y) can be expressed as follows: 

V(X <x,Y <y) = cxp \ - f ^ V ^/i(ds)| . (2) 

^ Js ^ y 
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The random vector in ([T]) is said to be standard, if it has standard 1-Frechet 
marginals or equivalently, Jg fdjj, = Jg gdfi — 1. It is weh known that any 
bivariate a-Frechet random vect or can be easily transfor med into a standard 
1-Frechet random vector (see e.g. IStoev and Taaau ( 20051) . Proposition 2.9). 



For the spectral representation of standard bivariate 1-Frechet distribu- 
tion (^j_S2ecifi£_choicesof(£j_2^ and {S, fj) often appear in the literature (see 



e.g. Ide Haan and Ferreiral (|2006[ ). Theorem 6.1.14). We will provide examples 
using the following one for our convenience. Set /(s) = 2s,g{s) = 2(1 — s) for 
s e S = [0,1], and — H, a probability distribution on [0, 1] with mean 1/2. 
In this case, 

P{X <x,Y <y)^exp(- — V ^^^—^H{ds)) . (3) 
^ Jq X y ) 

Many of our examples are constructed by choosing a specific H, and their joint 
cumulative distribution functions often do not have simple forms. 

The bivariate a-Frechet distributions arise as limits of i.i.d. bivariate reg- 
ularly varying random vectors, and it is ofte n convenien t to u se th e notion 
of vag ue convergence, denoted by 'A' (see e.g. IKallenberg ( 19861 ) and iResiiick 



(1987 1)). In particular, a random vector (F, VF) is said to be regularly varying 
with index a, if 

^F[(^,^)e.]^K-) (4) 

in Af+(E). Here, Kf e Wi/a, E = [0,oo]2 \ {(0,0)}, A/+(E) is the space of all 
nonnegative Radon measures on E, the limit i/ S E is non-zero and v{c-) = 
c~"i'(-), for all c > 0. The measure v is called the exponent measure, and it is 
said to have non-degenerate marginals, if J^((x, oo] x [0, oo]) and ^{[0, oo] x (y, oo 



are no n-degenerate in x and y, respectively. By Proposition 5.11 in iResnick 



([1987|), when a = 1 and v has non-degenerate normalized marginals, one can 



write 

2{l-w) 



/{([0,x]x[0,y])^}=^'^ 



y 



H{dw) 



such that H{dw) = J^2ujH{dw) — 1. In this case, letting {Vi,Wi),i = 
1, . . . , n be i.i.d. copies of {V, W), the vague convergence dH) is equivalent to 

(-yv.,-yw?i^{x,Y) (5) 

with {X,Y) defined as in ([2]). 

For any random vector {V,W), we say that V and W are asymptotically 
independent, if ^ holds with independent X and Y. This corresponds to the 
case when / and g have disjoint supports in ([T]), or H = {S{o} + (5{i})/2 in ([3]) 
(after normalization), or the exponent measure in (|4]) concentrates on the two 
axes {(x, 0), X > 0} U {{0,y),y > 0} with non-degenerate marginals. 
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Throughout this paper, for all measurable functions / and measurable sets 
So C S, we write \\f\\g^ = Jg^ \f{s)\fi{ds), and define 

Dt:={seS: f{s)/g{s) < t} and Et:^S\Dt. (6) 

Here and in the sequel, we will use the convention 1/0 = oo. For any standard 
bivariate 1-Frechet random vector {X,Y) as in ([1]), two important identities are 



lim II/IIb = / /l{g=o}d/i and lim 



\9\\d, = \\9\\s 



= 1 



which follow from the dominated convergence theorem. Furthermore, the joint 
distribution © can be expressed as 

<x,Y< y)^e^p[-[x-'\\f\\^^^^+y-%\\^^^^)] . (7) 

Finally, for any random variable Z and u > 0, we consider a thresholded version 
of Z denoted by 

Z{u) = u + {Z — = max(Z, u) (8) 

as m IZhand (|200^ ). We consider such thresholded random variables in order to 
analyze the quotient correlation coefficient ( (|2ip below), although it will turn 
out that the threshold value u does not play an essential role in the asymp- 
totic behavior (see Remark [5] below) . For a random vector {X,Y) in M^, let 
{{Xi,Yi)}i=i,,,,^n denote n i.i.d. copies of {X,Y). 



3 Ratios of Frechet Random Variables 



In this section, we provide an explicit formula for the distribution of X{u)/Y{u), 
u > 0. Here X{u) = X it u = 0. This leads to a necessary and sufficient 
condition for X{u)/Y{u) to be regularly varying. In the sequel, all the bivariate 
1-Frechet random vectors {X, Y) with representation ([1]) are assumed to be 
standard. We write at ~ bt if limt_>oo o-t/bt = 1- 

Our first result is an explicit formula for the joint distribution of (X/Y, Y). 



Proposition 1. Consider (X, Y) as in ([T]). Then, for alH > 0, u > 0, 



¥{X/Y <t,Y >u)= 1 



1 — cxp 



ll.9llz.,+<-'ll/ll 



In particular, when iJ.{Et) — 0, we have F{X/Y > t) 



(9) 



The proof borrows a result from IWeintraubl ( 199lh . and is deferred to Sec- 
tion [6l 

Remark 1. Prop osition [1] can be seen as a spe cial case of the conditional limit 
law established in Heffernan and Resnick ( 20071 ) (Propositions 4 and 5 therein), 
where (X, Y) are assumed to satisfy certain regular- variation type condition. 
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Their results describe the asymptotic hmit of t¥[{X /Y,Y /t) G •] as t — > oo. 
Here, thanks to the spectral representation, we provide an explicit formula for 
the joint distribution of {X/Y,Y). In particular, our result readily implies that 



lim P(X/r <t\Y >u) 



lim 



\X/Y <t,Y >u) 



— lim I 1 



p(y > u) 

11/11^^ ^-i l-exp[-(||g||^^+t-iH/H^Ju-i] 

't\\9\ 



Dt 



1 — exp(— u 1) 



\s\\d, = / 5l{/<t9}dM, 

JS 



which recovers Equation (32) in iHeffernan and Resnidd (|2007|). For 
a more general and geometric treatment of the conditional limit law, 



see 



Balkema and Embrechtd (|2007l) . 



Now, the distribution of X{u)/Y{u) follows as a corollary. 
Corollary 1. Consider {X,Y) as in Then, for &\\t>l,u> 0, 



\Y{u) 



>t = 1 



l/ll 



Et 



-iL / \\9\\D,+t-'\\f\\E, 
1 — cxp I 



In particular, for all u > 0, as t 

p( — > t 

\Y{u) 

Proof. Observe that for t > 1 



oo, 

t 



[1 - cxp(-^/-i)] 



(10) 



(11) 



.Y{u 

= p(x >tu,Y <u) + p(x > ty, y > u) 

= P(y <u)- P(X <tu,Y <u)+ p(y > u) - V{X <tY,Y > u) 
= 1 - V{X < tu, Y <u)- ¥{X <tY,Y >u). 

Plugging in ([Tj and ([9]) yields (flOl) . whence (fTTj) follows immediately, noting 
limt^,,||5||^^+i-i||/||^^ = = □ 



Inspired by pT|) . define 

7W := 



l/ll 



7 / /l{/>ts}dM- 



(12) 



Clearly, \imt^ooj{t) = 0. Then, provided j{t) > 0, or equivalently, fi{Et) > 
for all t £ (0,oo), (fTTj) implies that for all x > 0, 



P{X{u)/Y{u) > tx) 7(te)/(l + 7(te)) j{tx)_ 
P{X{u)/Y{u) > x) ^ 7(2:)/(1 + 7(.t)) ^ 7(2;) 

Therefore, we have thus proved the following result. 



as t 



00 . 
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Theorem 1. Consider {X,Y) given in ([1]) and suppose iJ.{Et) > for all t G 
(0,oo). Then, 



X{u) /Y{u) e RV-a for some (all) u > 0, if and only if 7 e RV-a ■ (13) 

In other words, for bivariate l-Frechet random vector {X,Y), to study the 
regular variation of X/Y, it is equivalent to study the regular variation of j{t) 
in p^ . based on the spectral functions. We will see in the next section that 
for many well-known examples, the regular variation of j{t) follows from simple 
calculations. 

Remark 2. When iJ.{Et) — 0, or equivalently ^{t) — Q for t large enough, we 
have F{X{u)/Y{u) > t) < F{X/Y > t) ^ hy Proposition [B This situation is 
relatively simple, and we do not study this case in this paper. 

Remark 3. A more general setting should be to consider (X, 1/Y ) in some do- 
main of attractions, such that X/Y £ RV-a- Maulik et al. ( 20021 ) investigated 
certain general cases under this framework. Our specific case, however, is not 
covered by their results. See more discussion in Section [T) 

From now on, we say {X, Y) has ratio tail index a, if (|13l) holds. We imme- 
diately have the following consequences. In the sequel, for any non-decreasing 
function U on R, let U^{y) := inf{s : U{s) > y} denote the left-continuous 
inverse of U. 

Corollary 2. Consider {X,Y) given in ([T]) and suppose fJ.{Et) > for all 
t € (0,00). Suppose that {X,Y) has ratio tail index a, then the following 
statements hold: 

(i) a — 1 if and only if ||/||^^ is slowly varying, 

(ii) for all u > and for all sequences k„ ^ (l/7)^(Cn), 

X,iu) 



1 

-V 



1 Y^{u) 



Co , (14) 



where C = 1— exp(— u ^) and Cq. is a standard a-Frechet random variables, 
and 

(iii) for any sequence {itnjngN such that u„ = o(n), the convergence in (jl4p 
holds with u replaced by u„ and k„ ^ (1/7)'*^ («/«„). 



Proof. Part (i) is trivial. Part (ii) follows from Proposition 1.11 in iResnick 
(|1987^ and part (iii) can be proved by a similar argument. □ 

An important consequence of Theorem [1] is that, the tail of the ratio is 
always lighter than or equivalent to the tails of X and Y. Indeed, (fT2|) implies 
7(i) — 0{t^^) as < — > 00, thus 7 e RV-a implies a > 1. The case a = 1 includes 
the case when X and Y are independent Frechet random variables. This follows 
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from iBreimaij (|1965[ ): if X £ RV-i, Z has finite 1 + e moment for some e > 0, 
and X and Z are independent, then XZ S RV-i (here, Z — 1/Y has standard 
exponential distribution) . 

However, for dependent X and Y, the ratio tail index a can still equal 1. 
By Part (i) of Corollary [21 a simple case of a = 1 is when limi^oo ll/l!^;^ = 
Jg,/l{g=o}dM > 0- This is the case, in view of the spectral representation 
when there is a point mass at {!}, or equivalently when one can write X = 
Xi V X2, such that Xi is 1-Frechet or equal to almost surely, X2 is 1-Frechet, 
and X2 is independent of Xi and Y. If in addition Xi = 0, then X and Y are 
independent. We can also have examples such that a = 1 and Jg fl^g^Qjdfi — 0, 
by constructing to be slowly varying, in a similar way as in Example [5] 

below. 

Next, we compare the ratio tail index and the tail dependence index defined 

by 

A = lim P{X >t\Y >t), (15) 



t- 



provided the limit exists. The random variables X and Y are asymptotically 
independent if A = 0, and asymptotically dependent if A € ( 0, 1]. The ta il 
dependence index has been widely studied and applied . See e.g. ISibuva ( 1960[ ). 



[de Haan and Resnick (1977h . iLedford and Tawii among others 



It is easy to see that when {X, Y) is bivariate 1-Frechet, X and Y are asymp- 
totically independent if and only if they are independent, as {X, Y) is in its own 
domain of attraction. As a consequence, the tail dependence index (when A > 0) 
measures the dependence strength of bivariate 1-Frechet random vectors. To 
understand the difference between the tail dependence index and the ratio tail 
index, observe that in our case, 

lim P(X > t |y > i) = |lmin(/,5)|L = / 2 min(s, 1 - .s)i7(ds) . 

The tail dependence index is determined by the spectral measure H on (0,1) 
(excluding {0, 1}), while the ratio tail index is determined by the behavior of H 
on a small neighborhood of {0, 1}. 

Example 1. Consider {X,Y) given by 

(X, Y) = ((1 - p)Z2 V pZs, pZi V (1 - p)Z2) , (16) 

where p e [0, 1] and Zi, Z2 and Z3 are independent standard 1-Frechet random 
variables. This corresponds, in ([3|), to choose 

^ = f ('^{0} + ^{\}) + (1 - P)'5{l/2} , 

where 6s^a} is the unit point mass at a. It is easy to see 7(i) = pjt for t > \. 
Therefore, XjY G that is, the ratio tail index always equals 1 for p > 0. 

Observe that in this model, however, the dependence strength between X 
and Y varies for different choices of p. Observe that P(X = Y) = P{Z2 > 
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p/(l — V Za)) = (1 — p)/(l + p), which increases from to 1 as p decreases 
from 1 to 0. By direct calculation, the tail dependence index equals X — 1 — p, 
which reflects the dependence strength. The ratio tail index does not. 

We conclude this section with a remark on role of the threshold value u. The 
ratio tail indices of several concrete bivariate distributions are calculated in the 
following section. 

Remark 4. Theorem[T]and Corollary[2]indicate that the threshold value u does 
not play an essential role in the asymptotic behavior of X{u)/Y{u). Indeed, 
from (jl4p . the limit, if it exists, is independent of u > 0, up to multiplicative 
constants. This suggests that letting u — cx) is more interesting as n — > oo, and 
part (iii) of Corollary [5] shows that this would change the convergence rate, but 
the limit is still the same. 

4 Examples 

In this section, we provide several examples on the ratio tail index. The first 
example is a concrete one where X/Y is not regularly varying. 

Example 2. We construct an example based on the spectral representation ([3]). 
We will find an H such that ^{t) ~ exp(— a.s t oo. Here exp(— t) is not 
regularly varying and therefore, by Theorem [1] X/Y is not regularly varying. 

Suppose H has Radon-Nikodym derivative h(s) w.r.t. the Lebesgue measure 
on [0, 1]. First, we need to find h{s) such that 



for t large enough and some constant C to be chosen later. Solving the above 
equation, we obtain 



One can choose a function /ig : [0, 1/2] K and tune C properly, such that 



yields a probability measure on [0, 1] with mean 1/2. (For a simple choice, one 
can consider H to be symmetric by setting ho{s) — hi{l — s) for s G [0, 1/2].) 
Thus, for {X, Y) corresponding to the spectral representation ([3]) with H given 
in (|17p . the extremes of X are compressed by the extremes of Y, resulting in an 
exponentially light tail of X/Y, which is not regularly varying. 

The next example is the spectrally discrete bivariate 1-Frechet random vec- 
tors. 






(17) 
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Example 3. A bivariate 1-Frechet random vector (X, Y) is spectrally discrete, 
if it has the foUowing representation 



(18) 



where Zi, . . . , Zm are i.i.d. standard 1-Frechet random variables, > 0, 6^ > 
0,1 = 1, . . . ,m. This model corresponds, in ([T]), to choose 5 = {1, . . . , m}, fi as 
the counting measure on S, f{i) = ai and g{i) = bi. 

Consider (jl8p with ai/bi strictly increasing. If am/bm < oo, then clearly we 
have linim^oo\/^=i Xi/Yi = am/bm almost surely. Otherwise, if am/bm = oo 
{bm = 0) then 7(t) = am/t for t > am-i/bm-i- Therefore, Theorem [1] implies 
X/Y G RV-i. That is, the ratio of spectrally discrete vectors is either bounded, 
or with tail index 1. 

Below, we calculate the ratio tail indices of two popular bivariate Frechet 
distributions, the logistic model and the mixed model. These models have 
been studied c areful ly and applied widely to many real data analyses. See 
Beirlant et al. I (12004 ' and the references therein for detailed results on these 



models and their modifications. 

Example 4 (The logistic model ( Gumbel ( 1960l ))). Consider the bivariate 1- 
Frechet distribution 



F{X <x,Y<y)^ exp 







1 






1 



,x > 0,y > 



(19) 



for a G [l,oo]. Here, the parameter a characterizes the dependence between X 
and Y. The random variables X and Y are independent, if a = 1, and they are 
fully dependent, i.e., F(X = Y) = 1 , if a = oo. 

Example 5.13 in iResnickl (|l987l ) showed that this distribution has a repre- 
sentation ([T|) with, when a > 1, 

f{s) = {a- l)."-i(l - s")-i/",g(.s) = {a- l)s«-2 , 

and {S,ii) = ([0, l],Leb). By straightforward calculation, we have 

Et^\se [0, l]:M>t]^\se [0, 1] : > t"] = ( 7 ^^T^, 



Therefore, writing 's — s 



f{s)ds 



(l + tQ)l/Q 



r ^(i-3)-v"d^a 



i V 1 + 



1-1/0 
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Theorem [T] implies that the logistic model ([T^ has ratio tail index a. This 
provides a probabilistic i nterpretation of a. Another way to interpret the pa- 
rameter a is provided bv iLedford and Tawn Therein, it is shown that 
1/a equals the limit probability that the component maxima do not occur at 
the same observation, i.e.. 



lim P(argmaXj^ 



argmaXj^ 



Remark 5 (Estimation of the logistic model). For the logist ic model in Ex- 
ample HI curiously one may want to apply the Hill estimator ( Hill ( 19751 )) on 
i.i.d. copies Ri '■— Xi/Yi to estimate a. Let R{i^n) < ■ ■ ■ < R{n,n) be the order 
statistics of {i?i},;=i....,r!- Then, the Hill estimator is defined by 



IH = jH{k,n) := 



1 ^ 

-^log i?(,,„) 



k 



log R(k 



+ l,n) I 



depending on a threshold integer value k. This estimator is consistent in the 
sense that 'Jh — » in probability, a s lon g as k/n — )■ and fc — >■ oo as 

n — >■ oo (see e.g. de Haan and Ferreira ( 20061 ) Theorem 3.2.2). Furthermore, 
the asyi nptotic normality of Jh is guaranteed by the second-order condition 
(see e.g. lde Haan and Ferreiral (j2006[) . Definition 2.3.1, Theorem 3.2.5). Indeed, 
since the explicit formula of F{X/Y > t) is available, after some calculation we 
can show 

Vk{jH - a-^) ^ Af{0,a-^) , 
for k = [n'^/2J , /3 G (0, 2/3). When ^ 2/3, the Hill estimator has the optimal 
rate (ri^/"^), but the limit will be a non-centered normal distribution. 

However, applying the Hill estimator here has little practical interest, since 
in this parametric model the maximum likelihood estimator works well with 
better rate {n^^"^): when a > 1, the estimation problem is regular and when 
a = 1, the n on-regu l ar be havior of the maximum likelihood estimator has been 
addressed bv lTawnl (Il988h . 

Example 5 (The mixed model ( Gumbel ( 19621 ))). Consider the bivariate 1- 
Frechet distribution 

k 



^(X <x,Y <y) ^ exp 



1 



y 



,x > 0,y > 



(20) 



for k G [0, 1]. Here, fc = corresponds to the independent case, though k ~ 1 
does not correspon d to the full d ependence case. By a similar calculation as 
in Example 5.13 in iResnickl ([1987|), this distribution has a spectral representa- 
tion ([3]) with 

i7(.) = fcLeb(.) + ^[<5{o}(-) + <5{i}(-)], 

where (5{o} and S^iy are unit point masses at {0} and {1}, respectively. Straight- 
forward calculation shows 



lit) 



l/ll 



1-k + k 



1 



t 



1 + t 
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Therefore, 

r (1 - /c)t-i if o<k<i 

When fc = 1, the ratio tail index equals 2, and 



i=l 



where C2 is standard 2-Frechet. When k £ (0, 1), the ratio tail index equals 1, 
and 



where Ci is standard 1-Frechet. The asymptotic behavior of the maxima of the 
ratios changes dramatically at /c = 1. 



5 The Quotient Correlation and Independence 
Test 

In this section, we apply our results on the ratio tail in dex to s tudy the asymp- 
totic behavior of a test statistic recently proposed by IZhand (|2Q08l ). Therein, 



Zhang proposed the gamma test for testing the independence between two ran- 
dom variables, based on their ratios. A similar test was also proposed aiming 
at testing the asymptotic independence (tail independence). The asymptotic 
behavior of the test statistics proposed have so far been studied only for in- 
dependent Frechet random variables. Here, we establish asymptotic results for 
jointly Frechet random variables with arbitrary dependence structure. We show 
that, in most of the dependence cases, the power of the hypothesis test goes to 
one as the samp le size increases to infinity. 

Zhang ( 2008[ ) essentially focused on the test statistics of the form 



Qu,n 



ma.-xi<i<„{Xi{u)/Yi{u)) + iiia-xi<i<n{Yi{u)/Xi{u)) - 2 
ma.xi<i<n{Xi{u)/Yi{u)) x maxi<i<n{Yi{u) / Xi{u)) - 1 ' 



where {Xi{u),Yi{u)) are i.i.d. copies of {X{u),Y{u)) (recaU ([8])). When u ~ 0, 
Qn '■= 90, n is called the quotient correlation coefficient and when u > 0, qu,n is 

called the tail quotient correlation co efficien t. 

When X and Y are independent, IZhand (l2Q08l) showed that 



Gvm^v|^)-((i---K!'',(i-e--Kr) 



and 



nqu,: 



(1 - e-V«)-i ( J_ + 1 ) A r(2, (1 _ e-i/")-i) , (22) 

Ci Ci 
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where C^^'' and C^^^ are independent standard 1-Frechet random variables. Here 
T(k,9) stands for the gamma distribution, which equals the distribution of the 
sum of k independent exponential random variables with mean 9. Recall that 
the inverse of a standard 1-Frechet random variable has standard exponential 
distribution. 

Based on_(f22|) with li = 0, a hypothesis testing for independence was designed 
m IZhand (| 2008[ ) with the following null and alternative hypothesis: 

Hq : X and Y are independent and Hi : X and Y are dependent. (23) 

The test statistic ng„ then has gamma limit distribution as (j22p under the null 
hypothesis, the test (|23p is thus referred to as the gamma test. When m > 0, a 
similar hypothesis test was designed for testing asymptotic independence. 

In this section, we address the asymptotic behavior of g„^„ for dependent 
bivariate 1-Frechet random vectors. Recall again that for jointly Frechet random 
variables, the independence and the asymptotic independence are equivalent. 
Therefore, we will focus on the independence test (|23|) and in particular g„ = 
go,n- This essentially requires to investigate the limit of the joint distributions 

1_ w X,{u) _1_ w Y,{u) 

with u = as n — >■ cx) for some suitable sequence {K^}„gN- However, all 
the asymptotic limits in this section would be the same for all w > 0, up to 
multiplicative constants depending on u. We choose u — also for the sake of 
simplicity. 

We first introduce some notations, since we need to deal with two groups 
of symbols, corresponding to X/Y and Y/X respectively. By default, the sym- 
bols with a sign '-|-' ('— ' resp.) correspond to the ratio X/Y (Y/X resp.). In 
particular, for {X,Y) as in ([1]), write -f+{t) = j{t) as in (fT2|) and 

^ /sgl{3>«/}dM 

7 U ^ 

We have shown in Section |4] how to calculate the regular variation of 7+(i) for 
several models. Here 7-(i) can be treated similarly. In particular, if /(s) = 
g{l — s) for s S [0, 1] and /i is symmetric on [0, 1], or equivalently iJ in ([3]) is 
symmetric on [0, 1], then 7+(i) = 7 - (t) . The follo wing theorem and its corollary 
generalize Theorems 3.1 and 5.3 in lZhangI ( 20081 ). 



Theorem 2. Consider {X,Y) given by Suppose n{Et) > and fi{Dt) > 
for all t e (0, oo). Then, j± £ RV-a± for some a+ > 0, a_ > 0, if and only if 

(^Vf'^Vi)-(C^.'C«J, (24) 



where ~ (l/7±)^(^) and (a± are independent standard a-t-Frechet random 
variables, respectively. When a+ = a- = a, (a+ and in ((24|) are interpreted 
as two independent standard a-Frechet random variables. 
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Proof. The 'if part follows from Tlieoreni[TJ We show the 'only if part. Write 
R+ = X/Y and R' = Y/X. Observe that, for > 1, 

P(i?+ < t+,R- < i_) = 1 - P(i?+ > <+) - V{R- > t_) . 

Therefore, for all t^^t^ > and n such that K^t^ > 1 and > 1, 

logP(^ V ^ < i+, V ^ < i_) ^ -n [P(i?+ > t+K+) + F{R- > 



? — 1 2—1 



By definition of k^, we have 

+ P(i?+ > 

and similarly nP(i?^ > t_K~) ^ Therefore 
n[F{R+ > t+K+)+P{R- > t^K-)] 



(25) 



+ P(i?+ > t+K+) 



P(i?+ > t+ + nF{R- > t+K-) t_ 



as n oo. Combined with p5|) we have thus proved (l24l) . □ 



As a corollary, we establish the asymptotic behavior of g„ . The test statistic 
qn has different asymptotic limits, depending on different regular-variation type 
behaviors of 7±. We write a„ = o(6„), if lim„_>.oo a„/6„ = 0. 

Corollary 3. Suppose that the assumptions in Theorem [2] hold. Consider 
4 (l/7±)^(n). Then, 

1 C 

K^g„ h - — if ^ Ck^ for some C > . (26) 

with Cq± as in Theorem [5] Note that a+ > a- implies = o(k~), and 
^ Ck,~ implies a+ = 



Proof. We only prove ((26| . The proofs for the other two cases are similar. For 
the sake of simplicity, write 

i?+ = V$' ^"=V# (27) 
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and i?± = Rt/i^t- Then, by dH]), we have 

Now, Theorem [5] imphes {R^,R^^) (Ca+ ; C'^^Ca- )■ By the continuous map- 
ping theorem, we have proved (pS)) . □ 

Corohary [3] provides theoretical support for the gamma test. Indeed, it 
shows that as long as ^ (l/7-|_)'*~(n) = o(n), then nq„ explodes quickly as 
n — >■ 00. This means in this case, the gamma test rejects the null hypothesis with 
probability going to one as the sample size increases to infinity. Furthermore, 
the following example shows that when ~ C±n for some constants C± > 0, 
the gamma test still performs reasonably well. In effect, this is indeed the 'worst' 
case that the gamma test could encounter, provided the ratio tail index exists. 

Example 6. Recall the model (fT6|) considered in Example [TJ 
{X, Y) = ((1 - p)Z2 V pZa, pZi V (1 - p)Z2) . 

In this case, 7±(i) = p/t for t > 1, and the tail dependence index A = 1 — p. 
Corollary |3] yields 



nq„ 



Now, observe that in the strong tail dependence case, i.e., p is close to 0, the 
power of the test is high. Consider a test of level 0.05. Let qp{/3),/3 £ [0,1] 
denote the lower /3 quantile of the distribution of r(2, p^^). For this model, the 
power converges to 

l-g;i(gi(0.95)). (29) 

Figure [1] illustrates the power of the test as a function of p € [0,1]. We see 
that the test performs reasonably well as long as the tail dependence index 
A = 1 — p is not too small. In addition, we also observe that the test statistic 
nqn converges quickly. For n = 20, the power of the test is already close to the 
limit one in 1^^. 

Remark 6. In this experiment, we used a modified gamma test instead of the 
original one. Namely, instead of g„ — qa^n defined as in (|2ip. we used a slightly 
different statistic qn R^R:^/{R^ + i?^) with R^ as in (P71) . We do not use 
the original test statistic qn because, in this example when p is small and the 
sample size n is small, with very high probability the go,n would equal 0/0, 
which is not well defined. 

One can show that nqn and nqn have the same asymptotic distribution. 
Indeed, if we replace q„ by g„ in Corollary [31 all the statements remain valid. 
The proof will be the same, except that we have K^gn = (i?^ + 
instead of 
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Power of the Gamma Test 




Figure 1: Power of the gamma test (solid line) in the 'worst scenario' (Exam- 
ple |6]) for selected p's with n = 20, based on 1000 simulations. The dashed line 
corresponds to the limit power calculated in (|29|) . 



Remark 7. The model (fT6| represents the 'worst' case that the gamma test 
can encounter. This is the case when the power of the gamma test is low, but 
the tail dependence is the strongest possible. 

To design such a scenario, we need rtg„ to converge and {X, Y) to have the 
largest tail dependence index possible. To have ng„ to converge, by Corollary[3] 
we must have ~ C±n, which is equivalent to the fact that H has point masses 
at {0} and {1}. Now, the tail dependence index is maximized by concentrating 
all the measure of H on (0, 1) at {1/2}. 

In this example, the gamma test performs poorly when p is close to 1. How- 
ever, one should not expect any independent test to perform well, as {X, Y) 
as in (|16p can be seen as a pair of independent random variables (pZi, pZ^), 
slightly perturbed by (1 — p)Z2 (via the 'max'(V) operation). 

InlZhangl (l2008h . other versions of qn were also proposed, aiming at dealing 



with arbitrary bivariate distributions. In principle, the data needs to be first 
transformed to have 1-Frechet marginals. It is an intriguing problem to char- 
acterize how the dependence structure changes during such a transformation. 
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The investigation along this hne is out of the scope of this paper. 

Remark 8. Observe that for continuous random variables X and Y with 
c.d.f. Fx and Fy, the joint distribution of (-1/ logFx(X), -1/ logi^y (F)) wiU 
have marginal standard 1-Frechet distributions. However, most of the time this 
transformation does not lead to a bivariate 1-Fre c het distribu tion. For studies 
on these distributions, see e.g. iLedford and Tawnl (|l99(l Il997h . 



6 Proof of Proposition [T] 



We borrow a result from IWeintraubl (Il99lh . Weintraub's work is based 
on the min-stable distribution (processes), which can be equivalently trans- 
formed into our setting. Namely, if {V, W) has a min- stable distr i bution with 
spec tral functions / and q g L^ iS,fi) according to IWeintraubl ( 1991 ) (see 
also Ide Haan and Picka^ (| 19861) ). then {X,Y) {1/V,1/W) is jointly 1- 
Frechet with represent ation in (ID). 

Now, Lemma 3.4 in IWeintraubl ( 199ll) becomes, for {X, Y) as in ([1]), 



P{X <x\Y = y) 



exp 



__5(f) 
X y 



r /(£i>-i. 

I 9(s) V I 
1 



y 



= ll.9llDt exp ( 
Then, letting Py denote the distribution of Y 

/>oo 

P{X/Y <t,Y >u) = / F(X <tr\Y- 

J U 

fix <tr\Y 



/^(ds 



3(s)l//(£l<-l.M(ds) 



(30) 



exp 



r)Py(dr) 



r ^ exp(— r ^)dr 

^ + t-'\\f\\E.-\V 



dr . 



1.911b, = 11.91 



We have thus obtained ([9]). The proof is 



Observe that 1 — 
complete. 

Remark 9. Careful re aders may find th e difference between our definition of 
Dt in ^ and the one in IWeintraubl ( 199ll) . Lemma 3.4 (where different symbols 
are used). Therein, Dt = {s € S : f{s)/g{s) < t}, with the '<' in ^ replaced 
by '<'. After reading the two paragraphs of the proof of Lemma 3.4, one should 
see that the c orrect defin i tion o f Dt is as in ^ with '<'. See also the proof of 
Lemma 3.5 in IWeintraubl (Il99ll) . where Lemma 3.4 was applied with '<'. 

An alternative way to see quickly it is not correct to choose the definition 
with '<' is given next. Note that, using '<' instead of '<' would only lead to a 
different formula when Uglljj^ (or ||/||^^, resp.) is discontinuous at some to > 0. 
Consider 

F{t) :=P{X/Y <t) ^P{X/Y <t,Y > 0) ^ ^* 



t\\9\ 



Dt 
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which has jumps at to such that n{s : f{s)/g{s) — to} > 0. In this case, the 
definition with '<' would cause F{t) to be right-discontinuous at to- Indeed, 
both 11/11^^ = Jg^ /d/i and \\g\\ij^ — Jjj^gdfi would be right-discontinuous at 
to- But F{t), as a cumulative distribution function, should be right-continuous, 
which is a contradiction. 



7 Discussion 



Recently, several extensions of the notion of regular variation have been in- 
troduced. The main motivation behind them is to study in more details the 
asymptotic behaviors th at are no t captu red by the standard multivariate regular 
variation. In particular. iResnickl (|2002l ) introduced the notion of hidden regular 
variation, characterizing the dependence structure when compone nts of the ran- 



dom v ectors are asymptotically independent but not independent: iMaulik et al 



( 2002h investigated the regular variation of products of random variables, and 
applied the result to model the network traffic. We discuss our results from 
these perspectives. 

Hidden regular variation. Recall the definition of regular variation in the lan- 
guage of vague convergence ([4]) . When the exponent measure v concentrates on 
the two axes, the two random variables V and W are asymptotically indepen- 
dent. In this case, J/ does not provide useful information to characterize different 
dependence structures for asymptotically independent random variables. More 
sophisticated models are needed. 

This problem was first investigated bv lLedford and Tawn ( 1996 . 1997 ). and 
their models later on were generalized under the framework of hidden regular 
variation by Resnick (2002, 2007, 2008). Hidden regular variation is present 
when the vague convergence dH) holds with i/ concentrated on the axes, i.e., V 
and W are asymptotically independent, and in addition there exists k° = o(Kt), 
such that 



tF 



in M+(Eo) with Eq = (0, oo] x (0,oo]. Intuitively, the notion of hidden regular 
variation involves normalizing the vector (V, W) by sequences of constants of 
smaller order (k° = o(«^t)) than required in Thus, ignoring the two axes, 
certain dependence structure might appear in the limit. This dependence struc- 
ture is not captured, therefore 'hidden', when the 'standard' rate is taken. 

In our case, when (X, Y) is a bivariate 1-Frechet random vector, we will see 
that {V, W) — {X/Y, Y) are asymptotically independent, and there is no hidden 
regular variation. We first look at the convergence in form of 



X/Y Y^ 

Kt ' t 



(31) 



in A/+(E) with Kt ^ {l/'y)^{n). Here, we are in a slightly different situation, 
as the normalizing sequences nt and t are of different rates most of the time. 
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In this case, we say X/Y and Y are asymptotically independent with different 
rates. The rates are the same, i.e., Kt ^ ct for some constant c S (0, oo), if and 
only if /l{g^o}dM > 0, by Corollary [21 part (i). Otherwise, Kt = o{t). 

We first show below that X/Y and Y are asymptotically independent (with 
different rates if {l/^)*~{t) — o{t)). Indeed, for r > 0, y > 0, Proposition [T] 
implies, writing t = Ktr for simplicity. 



:.t([0,r] X (y,oo]) = tp(^<r,^>y) 

\ Kf t / 



X/Y ^ Y 



1 — cxp I 



ty 



'■'^"^ + o(l) = i + o(l) , as < ^ oo . (32) 



y y 

It follows that for all rectangles E = [t*!, 2/i) x [r2, 2/2) in Eq, Ivcut^oo T^t{E) = 0. 
Recall also that by (HH) (with u = 0), t¥{X/Y/Kt > r) ~ r"". This together 
with (15^ yields the following result. 

Proposition 2. Let {X^Y) be as in ([1]) and suppose ((T3| holds. Then, for 
K( ~ (1/7)^ (t), the vague convergence ([31]) holds in M-|_(E) with concentrated 
on the two axes. In other words. 



1=1 ' i=l 

where c}^'^ and (^;[^'* are independent standard Frechet random variables with 
indices a and 1, respectively. 

Now, we examine the hidden regular variation of this model. Suppose that 
there exist sequences n° — o{{l / "f)^ {t)) and t° = o(i), such that 

.°(.):=ff[(^,^)e-l A.°(.) (33) 

in Af+(Eo) for some non-zero v° G A/+(Eo), i.e., {X/Y^Y) is hidden-regularly 
varying. Then, by a similar calculation as in I\i2\ . we have 

z.°((0,r]x(y,oo]) = ^||.9||^^„^(l + o(l)). 

^tV * 

We see that v° cannot converge in Eq. Therefore, we have proved the following. 
Proposition 3. There is no hidden regular variation for {X/Y^ Y). 



Regular variation of products. iMaulik et al.l ([20Q2[ ) proposed different general- 
izations of ([4]), for the purpose of characterizing the regular variation of the 
product of random variables. Their main results focus on two different situa- 
tions. 
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First, ifV e RV^a and 



{vc X G)(-) 



in M+(D),D = (OjCXi] x [0,oo], where Va_{{x,oi^ — x 
probability measure with G(0, oo) — 1. Then, 



(34) 



> and G is a 



{V, VW) 



in Af+(©) with v G M+(B), determined by and G (jMauhk et all (|2002l ). 
Theorem 2.1). This case can be seen as a generahzation of Breiman's theorem, 
since VW has an equivalent tail as V. Indeed, (p4|) is referred to as the (new) 
definition of asymptotic independence of V and W therein. 

Second, suppose V G RV-ay > ^ RV-aw some ay, a^/ > and V and 
W are asymptotically dependent with different rates, or equivalently. 



in M+(Eo) with i/((0,oo]^) > 0. Then, VW € In this 

situation, the multiplication with a random variable changes the tail behavior. 

Our asymptotic result on the bivariate Frechet random variables (Theorem[Tl 
with V = X and W = 1/^), does not fall into any of these situations. Indeed, 
if ([T3| holds with a > 1, then X/Y has a lighter tail than X, which differs from 
the first situation; on the other hand, 1/Y always has exponential distribution, 
which is not regularly varying, whence we are not in the second situation either. 
Therefore, the following question arises. 

Question 1. Let V and W be two nonnegative random variables. Suppose 
V G RV-a (but W may not be regularly varying). Provide a sufficient condition 
on V and W such that VW is regularly varying with tail index /3 > a > 0. 
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